Intermittency and structure functions in channel flow turbulence 
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We present a study of intermittency in a turbulent channel 
flow. Scaling exponents of longitudinal streamwise structure 
functions, Cp/C3i are used as quantitative indicators of inter- 
mittency. 

We find that, near the center of the channel the values of 
Cp/C3 up to p = 7 are consistent with the assumption of ho- 
mogeneous/isotropic turbulence. Moving towards the bound- 
aries, we observe a growth of intermittency which appears to 
be related to an intensified presence of ordered vortical struc- 
tures. In fact, the behaviour along the normal-to-wall direc- 
tion of suitably normalized scaling exponents shows a remark- 
able correlation with the local strength of the Reynolds stress 
and with the r.m.s. value of helicity density fluctuations. We 
argue that the clear transition in the nature of intermittency 
appearing in the region close to the wall, is related to a new 
length scale which becomes the relevant one for scaling in high 
shear flows. 



Spatio temporal intermittency is one of the most in- 
triguing properties of fluid dynamics turbulence. In- 
termittency can be described by means of the scaling 
behavior of the structure functions S p (r) built out of 
the velocity difference, namely S p (r) = (Sv(r) p ) where 
Sv(r) — 8v(r) ■ r/r. In the inertial range, S p (r) scale ho- 
mogeneously with exponents ( p , i.e. S p (r) oc r^ p . Inter- 
mittency reflects in anomalous values (i.e. different from 
the Kolmogorov 1941 prediction C, p — p/3) of the ( p (see 
Frisch §). Although many efforts have been devoted to 
the understanding of intermittency in homogeneous and 
isotropic turbulence, the case of wall turbulence domi- 
nated by very strong shear flow is still under debate. In 
fact, while the decrease of scaling exponents towards the 
wall has been recently pointed out by both experimental 
and numerical investigations j^Jf , a physical explanation 
of this effect is still missing. 

The main objective of this paper is to analyse the scal- 
ing exponents (if any) and to explore their functional de- 
pendence on the non-homogeneous coordinate. To this 
purpose, we investigate numerically a channel flow 
possibly the simplest instance of a shear-dominated flow. 

The problem of characterizing the complex phenom- 



ena arising in the near-wall region, and their relation to 
the lack of isotropy has been analyzed in depth by sev- 
eral authors (e.g., see Antonia et al. [|| and L'vov and 
Procaccia Qj). With specific reference to intermittency, 
Kutznetsov et al. Js) presented an experimental investi- 
gation of fine-scale structure of turbulence for different 
shear-flows. Even though they recognize that scaling ex- 
ponents may be different for various flows or various lo- 
cations of the same flow, they do not seem to address 
the issue of the spatial dependence of the exponents £ p . 
This issue is examined in a recent paper by Camussi et 
al. jo) , in which the authors propose a technique to iden- 
tify coherent structures and relate them to the increase 
of intermittency in (homogeneus and non-homogeneus) 
grid turbulence. The present work is fairly distinct in 
purpose, since we analyse the spatial behaviour of the 
scaling exponents and their link to coarse-grained fea- 
tures of the flow in the near-wall region. Besides the 
fundamental interest on its own, the existence of such 
correlation could prove very valuable for the design of 
more efficent Large-Eddy-Simulation models (see for in- 
stance Scotti and Menevau JlO||). 

We have performed a direct numerical simulation 
achieving a high statistical accuracy (about 10 3 in time 
units Uo/h, where Uq is the centerline velocity and h 
is the channel half- width). Numerical simulations have 
been performed on a massively parallel machine using 
a LBE (Lattice Boltzmann Equation) code. The spa- 
tial resolution of the simulation was 256 x 128 x 128 grid 
points. Periodic boundary conditions were imposed along 
the streamwise (x) and spanwise (z) directions, whereas 
no-slip boundary conditions were applied at the top and 
the bottom planes (normal to wall direction, y). The 
Reynolds number is Re ~ 3000. Further details about 
the numerical scheme can be found in [jxijl and references 
therein. In the following we use wall units defined as 
y + = y ■ v* jv and v + — v/v* where v* is the friction 
velocity fl^ ]. In these units, the channel is 640 long, 320 
wide and 320 high. 

To study intermittency in the channel, we introduce 
the following y-dependent longitudinal streamwise struc- 
ture functions: 
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S P (r + ,y + ) = (\(v x (x + + r + ,y + ,z + ) - v x {x + ,y + , z + )\ p ) 

(1) 

The average is taken at a fixed y + value (the normal 
to wall coordinate). The quantities S p (r + ,y + ) have been 
measured for each value of y. Our data set allows enough 
statistical accuracy to estimate S p (r + , y + ) forp < 7. Due 
to the low Reynolds number, we use ESS [H| in order to 
extract C, p values. We remind that ESS consists of mea- 
suring structure functions as a function of S3 rather than 
in terms of space separation r. This procedure allows a 
much better accuracy for the evaluation of the scaling 
exponents, although it does not provide any estimate of 
the y + -dependence of £3- In order to compute the scaling 
exponents ( p (y) we have analyzed the ESS local slopes 
D M (r+,y+) = d\og(S p (r+,y+))/d\og(S q (r+,y+)) for 
each value of the y + coordinate. We have found two re- 
gions in y + , hereafter referred to as region H ("Homoge- 
neous") and region B ("Boundary") respectively, where 
well defined constant local slopes for the scaling expo- 
nents can be detected. Region H is close to the center 
of the channel (y + > 100) while region B is close to the 
viscous sublayer (20 <y + < 50). In region H, the scaling 
exponents ( P {H) are found to be approximately the same 
as the ones measured in homogeneous and isotropic tur- 
bulence. On the other hand, in region B the scaling ex- 
ponents ( P (B) have been found to be much smaller than 
C, P (H). Moreover, while in region H the scaling range 
starts at r + > 25, in region B the scaling range starts at 
r + > 50, consistently with previous findings H]. In the 
intermediate region between region H and region B, it is 
difficult to identify a range in r where a scaling exponent 
can be defined with enough confidence. 

In order to clarify the discussion, we show in Fig. ^ 
and Fig. ||the local slopes D$^{r + , y + ) and D^2{r + , y + ) 
respectively for y + = 30, 70, 80, 150. In the intermediate 
region (i.e. y + = 70,80), the analysis in terms of local 
slope does not provide a well defined scaling exponent 
since the plateau in r + is very short. In Fig. |^ and Fig. 
H we show C6/C3 an d C4/C2 respectively as a function of 
y + with the associated error bars. The large error bars 
in the region 50 < y + < 100 indicate that scaling expo- 
nents defined through D p ^ q (r, y + ) are poorly defined and 
should be considered just as effective exponents obtained 
by the power law fit of the ESS analysis. The situation 
described in Fig. |l]and Fig. || is similar for all the scaling 
exponents £ p computed in our analysis. Finally, in Table 
|| we list the numerical values of the scaling exponents 
for region B, region H and for the homogenous isotropic 
turbulence [|l5| . Our results indicate that there is a tran- 
sition in the nature of intermittency between region H 
and region B. While in region H intermittency is close 
to what has been observed in isotropic and homogenous 
turbulence, in region B much stronger intermittency is 
observed, which reflects in lower values of the scaling ex- 



ponents £ p for p > 3 and larger ones for p < 3. Inbetween 
the two regions, a competition between the two types of 
intermittency should take place, leading to a poorly de- 
fined scaling law. 

An important question to be addressed concerns the 
physical mechanisms which produce much stronger in- 
termittency in region B with respect to region H. A pre- 
liminary answer to this question is given in Fig. [5] where 
the momentum flux (v' x v'y) (a particular component of 
the Reynolds stress tensor) normalized by its maximum 
(v' x v'y) , is shown as a function of y + . The momentum 
flux has a peak within the region B while it goes to zero 
in region H. This behaviour clearly indicates that inter- 
mittency grows dramatically in region characterized by 
strong mean shear (strong momentum fluxes). The link 
between momentum flux and the increase of intermit- 
tency can be further investigated by looking at the inset 
of Fig. H where we represent Ce(?/)/C3(y) from Fig. || and 
the rescaled expression 

( e (y + ) _ CeW = ( Ce(B) _ te(H) \ (v' x v' y ) m 
Cs(y+) UH) \Cs(B) ( 3 (H)J (v' x v' y ) M 1 ' 

From Fig. ^, we can argue that the increase of inter- 
mittency should be related to the increase of momentum 
flux and, therefore, to the mean (local) shear. 
Moreover, it is well known that turbulent flows near the 
wall are characterized by well defined coherent struc- 
tures. In Fig. [5] we show the r.m.s. helicity h (h rms = 
((tl> • v)) rms ) as a function of y + which is again peaked in 
region B. Coherent structures carry a significant amount 
of helicity while dissipation is found to peak in the in- 
terstitial region between helicity-carrying structures |l6| . 
Indeed a clearcut anticorrelation between helicity fluctu- 
ations and dissipation is systematically detected in our 
numerical simulations. Thus, the alternate presence of 
regions of high helicity and regions of high dissipation 
may be responsible for the enhancement of intermittency 
in region B. 

A more quantitative way to investigate the increase of 
intermittency in region B, can be achieved by the fol- 
lowing argument. According to Howarth-Von Karman- 
Kolmogorov equation for homogeneous shear flows tur- 
bulence (see Hinze JlTj]), one can define a length scale 
L s (y) in terms of the mean energy dissipation e(y) and 
the mean shear E(y), as follows: 

In presence of mean shear E, for any scale r we can define 
two characteristic timescales, namely the eddy turnover 
time r/Sv(r) and 1/E. We expect that when the mean 
shear is large enough, the eddy turnover time is not the 
relevant timescale for energy transfer from large to small 
scales. The inequality r/5v(r) < 1/E gives the range of 
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scales r where the effect of shear should not be relevant to 
small scale statistics. By using the Kolmogorov estimate 
Sv(r) ~ e 1 / 3 r 1 / 3 , we find that the above inequality can 
be written as r < L s . Thus, for r < L s one expects that 
the scaling properties of turbulence are not affected by 
the mean shear. On the other hand, for r > L s one ex- 
pects that the mean shear may significantly change the 
amount of intermittency. In our numerical simulation 
L 8 {y) becomes small only in the region B, i.e. in the re- 
gion where an increase of intermittency is observed. This 
result confirms our finding of a rather clear transition in 
the physical nature of intermittency, somehow similar to 
the Bolgiano scaling appearing in the thermal turbulence. 
Hence, this preliminary analysis seems to indicate that 
the change of scaling exponents cannot be reduced to a 
perturbative effect in terms of the mean shear. As a fi- 
nal observation, by inspecting the peak value approached 
by the ( p {y) exponents near the wall, we find an interest- 
ing similarity with the values, £p S , pertaining to a passive 
scalar advected by a turbulent 3D homogeneous/isotropic 
velocity field [Q . Values of passive scalar s are shown 
together with our present data in Table |. 

Table | suggests that the passive scalar behaviour can 
be traced to these helicity-carrying coherent structures 
being passively advected by the flow. This observation 
is in qualitative agreement with the results reported by 
Pumir and Shraiman Jl9| . 

In conclusion, we have rather clear evidence that, in 
wall bounded turbulence, the increase of intermittency 
near the wall is strongly related to the increase of the 
mean shear. We have introduced a characteristic length 
scale L s induced by the mean shear whose physical mean- 
ing is equivalent to the Bolgiano scale for natural con- 
vection. Velocity fluctuations at scales r > L s are ob- 
served to be more intermittent than in homogeneous and 
isotropic turbulence. 
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TABLE I. Normalized scaling exponents in regions H and 
B for the present simulation, as compared with the values 
for homogeneous isotropic turbulence Cp° m and for a passive 
scalar Q, B . 



p 


(Cp/Cs) 1 "™ 


(Cp/C 3 ) h 


(Cp/Cs) b 


Sp 


Cp P 7C 3 PS 


1 


0.36 


0.37 


0.44 


0.37 


0.46 


2 


0.70 


0.70 


0.77 


0.62 


0.77 


3 


1.00 


1.00 


1.00 


0.80 


1.00 


4 


1.28 


1.28 


1.17 


0.94 


1.17 


5 


1.54 


1.54 


1.31 


1.04 


1.30 


6 


1.78 


1.78 


1.44 


1.12 


1.40 


7 


2.00 


2.00 


1.55 


1.20 


1.50 
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FIG. 1. Local slope D6,3 for four different values of y + . 
The straight lines represent the fit for H and B regions. 



FIG. 2. Local slope 1)4,2 for four different values of y 
The straight lines represent the fit for H and B regions. 



FIG. 3. Scaling exponents C6/C3 as a function of y + . The 
straight lines represent the fit for H and B regions. 



FIG. 4. Scaling exponents C4/C2 as a function of y + . The 
straight lines represent the fit for H and B regions. 



FIG. 5. The momentum flux (v x v y ) (dashed line), normal- 
ized by its maximum, is presented together with normalized 
dissipation (straight line) and r.m.s. helicity (dotted line). In 
the inset the scaling exponent £e / C3 i s presented together with 
the rescaled momentum flux as defined by equation 2 (con- 
tinuous line). 



FIG. 6. The characteristic scale L s (y + ) is presented as a 
function of y + (continuous line). The dotted line is reported 
to enlight the growth of L s out of the wall region. 
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